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Introduction 


The  primary  objective  of  this  work  is  to  develop  new  techniques  in  single  photon  emission  computed 
tomography  (SPECT)  breast  cancer  imaging  providing  improved  detection  and  characterization  of 
early  stage  breast  cancer.  Conventional  SPECT  breast  imaging  uses  a  horizontal  axis  of  rotation  and 
thus  requires  a  large  radius  of  rotation.  By  using  a  vertical  axis  of  rotation  (VAOR)  and  various 
camera  head  and  collimator  configurations,  the  effective  radius  of  rotation  (ROR)  can  be  reduced 
and  background  myocardial  and  liver  interference  can  be  minimized.  A  tiltable  head  SPECT  (TH- 
SPECT)  system  with  a  VAOR  further  allows  for  improved  imaging  of  lesions  near  the  chest  wall 
and  reduction  of  effective  ROR  [2].  Investigation  of  new  reconstruction  algorithms  for  the  TH- 
SPECT  geometry  is  then  needed  to  accurately  (and  efficiently)  reconstruct  TH-SPECT  data  [1]. 
This  research  investigates  the  feasibility  of  TH-SPECT  acquisition  and  reconstruction. 

Body  (outlined  with  Statement  of  Work) 

Task  1.  Develop  TH-SPECT  image  reconstruction  algorithm 

Refer  to  attached: 

[1]  BC  Pieper,  JE  Bowsher,  MP  Tornai,  J  Peter,  CN  Archer,  RJ  Jaszczak.  "Parallel-beam  Tilted-head 
Analytic  SPECT  Reconstruction."  IEEE  Transactions  on  Nuclear  Science,  in  press,  2002 


Task  2.  Investigate  various  VAOR  system  configurations  as  well  has  different  camera  head 

geometries  and  collimator  configurations. 

Refer  to  attached: 

[2]  BC  Pieper,  JE  Bowsher,  MP  Tornai,  J  Peter,  RJ  Jaszczak.  "Breast  Tumor  Imaging  Using  Tiltable 
Head  SPECT  Cameras."  IEEE  Transactions  on  Nuclear  Science,  vol.  48,  no.  4.  pp  1477-1482, 2001. 


Key  Research  Accomplishments 

■  Development  and  evaluation  of  FBP  TH-SPECT  image  reconstruction  algorithm  [1]. 

■  Comparison  of  FBP  TH-SPECT  algorithm  with  traditional  FBP  and  OSEM 
reconstruction  methods  [1]. 

■  Acquisition,  evaluation,  and  comparison  of  TH-SPECT  data,  including  contrast  and 
S/N  analysis  using  experimental  phantoms  [2]. 
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Reportable  Outcomes 


[1]  BC  Pieper,  JE  Bowsher,  MP  Tornai,  J  Peter,  CN  Archer,  RJ  Jaszczak.  "Parallel-beam  Tilted-head  Analytic 
SPECT  Reconstruction."  IEEE  Transactions  on  Nuclear  Science,  in  press,  2002 

[2]  BC  Pieper,  JE  Bowsher,  MP  Tornai,  J  Peter,  RJ  Jaszczak.  "Breast  Tumor  Imaging  Using  Tiltable  Head 
SPECT  Cameras."  IEEE  Transactions  on  Nuclear  Science,  vol.  48,  no.  4.  pp  1477-1482, 2001. 

Brett  Pieper,  Master  of  Science,  Biomedical  Engineering,  Duke  University,  May  2002 

Conclusions 

TH-SPECT  is  an  exciting  new  technique  in  SPECT  breast  cancer  imaging.  Focus  on  TH-SPECT 
feasibility,  especially  with  small  field-of-view  SPECT  cameras,  as  well  as  improved  reconstruction 
methods  may  someday  supplant  traditional  breast  cancer  detection  methods.  Initial  results  of  this 
research  have  indicated  the  advantages  of  TH-SPECT  and  successful  TH-SPECT  (analytic  and 
iterative)  reconstruction.  Further  research,  particularly  patient  studies,  will  be  required  to  realize  the 
full  potential  of  TH-SPECT  breast  cancer  imaging. 
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Parallel-beam  Tilted-head  Analytic  SPECT 
Reconstruction:  Derivation  and  Comparison  with 

OSEM 

Brett  C.  Pieper,  Member,  IEEE,  James  E.  Bowsher,  Member,  IEEE,  Martin  P.  Tomai,  Member,  IEEE,  Caryl 
N.  Archer,  Member,  IEEE,  Ronald  J.  Jaszczak,  Fellow,  IEEE 


Abstract —  Parallel-beam  Tilted-Head  SPECT  (TH-SPECT)  was  previ¬ 
ously  implemented  on  a  SPECT  system  for  its  potential  to  image  breast 
lesions  and  nearby  axilla  of  seated,  upright  women.  All  TH-SPECT  re¬ 
constructions  will  contain  artifacts  since  the  tilted  orbit  does  not  satisfy 
the  Orlov  sampling  criteria.  However,  it  is  not  dear  which  reconstruc¬ 
tion  method,  if  any,  is  better  suited  for  TH-SPECT  data.  Here  a  geomet¬ 
ric  derivation  of  the  ramp  filter  for  tilted  parallel-beam  geometries  is  pre¬ 
sented.  A  filtered  backprojection  (FBP)  algorithm,  using  this  filter,  was 
then  implemented  and  compared  with  an  iterative  Ordered  Subsets  Expec¬ 
tation  Maximization  (OSEM)  algorithm,  using  TH-SPECT  data.  A  breast 
scan  at  various  tilt  angles  was  simulated  and  used  to  generate  a  noise  vs. 
bias  study  for  both  methods.  Contrast  and  SNR  values  as  well  as  axial  elon¬ 
gation  present  in  all  TH-SPECT  reconstructions  were  characterized  using 
a  mini-Defrise  disk  phantom  placed  inside  a  tillable  breast  phantom  and 
imaged  from  0-15  degree  head  tilt.  A  tillable  breast  phantom  containing 
lesions  was  also  imaged  with  a  system  dedicated  to  prone  breast  SPECT 
from  0-30  degrees  to  evaluate  the  effects  of  incomplete  sampling  due  to 
greater  tilt  angles.  FBP  noise  vs.  bias  studies  indicated  a  greater  increase 
in  bias  with  tilt  angle  compared  to  OSEM  reconstructions.  At  small  tilt  an¬ 
gles  about  the  mini-Defrise  disk  phantom,  poorer  contrasts  were  obtained 
with  FBP  compared  to  OSEM  at  similar  noise  levels.  All  reconstructions 
of  the  tillable  breast  phantom  indicated  axial  elongation  at  greater  tilt  an¬ 
gles,  although  FBP  reconstructions  displayed  an  increase  in  stretching  dis¬ 
tortions  of  the  breast.  OSEM  SNR  and  contrast  values  were  higher  at  all 
degrees  of  tilt.  In  conclusion,  measured  results  indicate  OSEM  TH-SPECT 
reconstruction  provides  better  contrast  and  SNR  values  and  may  offer  bet¬ 
ter  shape  and  uniform  activity  distribution  of  the  breast  compared  to  FBP 
methods. 


I.  INTRODUCTION 

Tilted-Head  Single  Photon  Emission  Computed  Tomography 
(TH-SPECT)  parallel-beam  imaging  of  breast  lesions  has  been 
shown  to  offer  superior  lesion  contrast  and  signal-to-noise  ratio 
(SNR)  values  compared  to  conventional  SPECT  geometries  [1]. 
Reconstruction  of  TH-SPECT  data  requires  a  three-dimensional 
reconstruction  algorithm.  Iterative  methods  have  thus  far  been 
used  to  accurately  model  the  complex  (albeit  incomplete)  sam¬ 
pling  of  frequency  space  associated  with  the  tilted  circular  orbit. 
It  may  be  useful,  however,  to  reconstruct  TH-SPECT  data  using 
analytic  filtered  backprojection  (FBP)  methods  given  the  speed 
and  linearity  associated  with  such  methods. 

Manuscript  received  November  5,  2001.  This  work  was  funded  by  grant 
DAMD 17-01 -1-0226  from  the  Department  of  Defense,  PHS  grants  R01- 
CA76006  and  R01-CA33541  from  the  National  Institutes  of  Health,  grant  RG- 
99-0305  from  the  Whitaker  Foundation,  and  grant  DE-FG02-96ER62150  from 
the  Department  of  Energy. 

B.C.  Pieper,  M.  P.  Tomai,  C.N.  Archer,  and  R.  J.  Jaszczak  are  with  the  De¬ 
partment  of  Biomedical  Engineering,  Duke  University  and  the  Department  of 
Radiology,  Duke  University  Medical  Center.  J.  E.  Bowsher  is  with  the  Depart¬ 
ment  of  Radiology,  Duke  University  Medical  Center,  Durham,  NC  27710  USA 
(telephone:  919-684-7687,  919-684-7788,  919-684-7791,  919-684-7766,  919- 
684-7785,  e-mail:  pieper@dec3.mc.duke.edu,  jeb@dec3.mc.duke.edu,  mar- 
tin.tomai@duke.edu,  cnb3@duke.edu,  rjj@dec3.mc.duke.edu). 


Conventional  2-D  FBP  SPECT  reconstruction  for  parallel- 
beam  collimation  relies  on  the  standard  one-dimensional  ramp 
filter.  This  ramp  filter  has  also  been  shown  to  be  the  backpro¬ 
jection  filter  for  3-D  tilted  parallel-beam  geometries  [2].  Due 
to  insufficient  Orlov  Sphere  sampling  [3],  TH-SPECT  acqui¬ 
sitions  introduce  a  null  space  of  unmeasured  frequency  com¬ 
ponents  whose  frequency-space  volume  is  dependent  upon  tilt 
angle.  These  unmeasured  frequency  components  may  have  neg¬ 
ative  effects  on  image  reconstruction.  One  such  effect  is  axial 
stretching  which  produces  unwanted  distortion  of  the  patient’s 
body  contour  [1].  Another  related  effect  is  inaccurate  activ¬ 
ity  distribution  specifically  due  to  insufficient  sampling  of  the 
Orlov  Sphere.  With  these  negative  effects  in  mind,  it  is  un¬ 
clear  which  image  reconstruction  method  is  better  suited  for 
TH-SPECT  data.  This  paper  develops  a  geometric  derivation 
of  the  ramp  filter  for  TH-SPECT  and  then  compares  parallel- 
beam  TH-SPECT  FBP  reconstruction  with  an  existing  iterative 
Ordered  Subsets  Expectation  Maximization  (OSEM)  algorithm 
[1],  [4].  Simulation  studies  as  well  as  phantom  studies  are  per¬ 
formed  and  evaluated. 

II.  Methods 

A.  TH-SPECT  reconstruction 

Three  spatial  domain  coordinate  systems  are  used  (Fig.  1) 
and  are  defined  as  (x,y,z),  (x',y/,z!)i  and  (x",/'3z")-  The  ob¬ 
ject  is  expressed  within  die  fixed  (x,y,z)  coordinate  system  as 
f(x,y,z).  The  (x \yr,z!)  system  is  a  rotational  coordinate  sys¬ 
tem  rotated  about  the  z  axis  by  angle  0.  The  third  coordinate 
system  (x "j/'jz")  represents  the  tilt  of  the  system 

about  the  /  axis  by  tilt  angle  <J>.  The  object’s  frequency  domain 
representation  is  expressed  within  the  fixed  (/z,v,  t)  system  as 
F(/i,  v,x),  and  rotations  within  this  system  are  similarly  defined 
as  and  (/i",v", t"). 

Rotation  of  (x ,y,z),  or  similarly  (/a,  v,t),  by  angle  0  produces 
the  following  transformations 

y!  =  xcos0  +  ysin0 

y  —  -jcsin0  +  ycos0  (1) 

z  —  z. 

Likewise,  rotation  of  (x/,y',z!),  or  similarly  (j a'jV'jT'),  by  an¬ 
gle  <|>  produces  the  following  transformations 
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Fig.  1.  Illustration  of  tilted  geometry  and  (x,y,z)  coordinate  system.  Lower- 
right  illustration  indicates  breast  positioning  within  the  tilted  geometry, 
where  the  z  axis  passes  through  the  nipple  and  the  center  of  the  base  of 
the  breast.  Dashed  lines  indicate  axis  is  behind  object. 


jJf  =  xcos0cos<|>-Fysin0cos<t>  —  zsintj) 

yn  =  —  xsin0  +  ycos0  (2) 

z  =  Jccos0sin(l)  +  ysin0sin(()-hzcos(|). 


A.  1  Backprojection  filter  for  non-tilted  parallel-beam  geometry 

In  order  to  suggest  a  filter  for  TH-SPECT  data,  we  first  ex¬ 
amine  filtering  for  non-tilted  parallel-beam  geometries.  The 
Fourier  projection-slice  theorem  implies  that  by  taking  projec¬ 
tions  at  multiple  angles  0  with  <j>  =  0,  values  of  F(/z,v, x)  will 
be  known  along  radial  planes  centered  at  the  origin  and  perpen- 
dicular  to  the  projections  [5],  [6],  [7],  [8],  [9],  [10].  Figure  2 
(top)  illustrates  the  frequency  domain  representation  of  a  sin¬ 
gle  non-tilted  projection  at  angle  0,  with  co  being  the  distance 
from  one  arbitrary  sample  point  to  the  axis  of  rotation  (AOR)  t. 
Rotation  about  x  fills  F(fi,v,x)  with  multiple  circles  of  discrete 
sample  points.  Figure  2  (bottom)  illustrates  a  slice  at  a  constant 
value  of  x.  As  the  radial  distance  co  increases,  the  density  of 
sample  points  within  a  given  rotation  circle  decreases.  Specifi¬ 
cally,  the  density  of  sample  points  within  each  circle  is  directly 
proportional  to  the  inverse  of  the  circles'  circumference,  or, 
in  our  rotated  coordinate  system,  ^7.  Because  the  distance,  y, 
between  all  rotation  circles  is  constant,  the  density  of  circles 
within  the  plane  is  also  constant.  The  density  of  sample  points, 
which  can  be  expressed  as  the  density  of  circles  times  the  den¬ 
sity  of  points  per  circle,  is  therefore  proportional  to  The  fil¬ 
ter  co,  which  is  the  well-known  ramp  filter,  compensates  for  this 
varying  sampling  density.  Because  lower  frequencies  are  sam¬ 
pled  more  densely  than  higher  frequencies,  a  filter  is  needed  to 
correct  for  this  difference  in  sampling.  Because  each  slice  along 
x  is  identical  in  sample  density,  the  backprojection  filter  for  non- 
tilted  parallel-beam  data  is  one-dimensional  across  projections. 

A. 2  Backprojection  filter  for  tilted  parallel-beam  geometry 

For  tilted  parallel-beam  data  (<()  /  0),  Figure  3  (top)  illustrates 
the  region  of  frequency  space  sampled  by  a  projection  at  angles 


Axis  of  rotation 


Fig.  2.  (Top)  Illustration  of  Fourier  projection-slice  theorem  for  non-tilted  data. 
The  plane  indicates  the  frequencies  sampled  by  a  projection  at  angle  0.  to 
represents  the  distance  between  an  arbitrary  sample  point  and  the  x  axis. 
(Bottom)  Frequency  domain  slice  at  a  constant  value  of  x  showing  the  dis¬ 
crete  sample  points  at  distances  v"  =  v*  =  to  from  the  AOR  and  rotated  by 
0.  Distance  y  represents  the  constant  distance  between  circles  of  discrete 
sample  points. 


0  and  <)).  The  distance  0)  from  the  axis  of  rotation  x  to  a  sample 
point  is  0)  =  \J v//2  +  x"2  sin2  (|). 

Sampled  frequencies  correspond  to  discrete  values  of  v",  and 
each  value  of  v"  lies  on  a  circle  of  sample  points  having  a  ra¬ 
dius  (0,  with  one  sample  point  per  angle  0.  The  density  of  points 
along  a  given  circle  is  inversely  proportional  to  the  circles'  ra¬ 
dius,  that  is,  proportional  to 


yj v//2  4-  x//2  sin2 


(3) 


In  contrast  to  non-tilted  data,  the  distance,  y,  between  circles 
in  this  plane  is  not  constant  but  rather  depends  upon  v"  and  x" 
(Fig.  3,  bottom).  The  density  of  circles  is  proportional  to  the 
reciprocal  of  the  distance  between  circles,  that  is,  proportional 
to 
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for  a  projection 


0  =  90° 


Fig.  3.  (Top)  Illustration  of  Fourier  projection-slice  theorem  for  tilted  data. 
Plane  illustrates  the  frequency  space  sampling  for  a  projection  at  angles  0 
and  4>-  <*)  represents  distance  between  arbitrary  sample  point  and  the  x  axis. 
(Bottom)  Frequency  domain  slice  at  constant  value  x  ^  0  built  up  from  dis¬ 
crete  sample  points  along  v"  and  over  rotation  angles  0.  Frequencies  are 
not  measured  within  the  concentric,  inner  circle  of  radius  of  xsin<|>.  Un¬ 
measured  frequency  region  decreases  as  x  — >  0  and  becomes  Figure  2  (bot¬ 
tom)  when  x  =  0.  Note  distance  co  is  now  a  function  of  both  v"  and  xsinq. 
Distance  y  represents  the  distance  between  circles  built  up  from  discrete 
sample  points.  Note  circles  are  more  dense  near  o>  =  xsin<t>  and  less  dense 
with  increasing  co.  For  simplicity,  sample  points  were  only  shown  at  45° 
increments.  Conventional,  finer  sampling  would  diminish  this  non-uniform 
density  appearance  within  each  circle. 


As  v"  increases,  the  distance  between  neighboring  circles  in- 
creases.  These  circles  are  more  dense  near  to  =  x"  sin(|>  than  at 
higher  co.  For  co  <  |x"  sincj)|,  the  density  of  sample  points  within 
the  plane  is  zero.  For  w  >  |x"  sin<|)|,  the  density  of  sample  points 
within  the  plane  is  the  product  of  (3)  and  (4),  which  is  ^y. 

Finally,  because  the  spacing  of  sample  points  in  the  x  direc¬ 
tion  for  tilted  data  is  constant,  although  compressed,  a  third  con¬ 
stant  component,  can  be  included. 

Multiplying  these  three  density  components  yields 


Fig.  4.  Illustration  of  missing  conical  regions  present  in  all  TH-SPECT  fre¬ 
quency  space  representations.  Volume  of  missing  region  increases  with  tilt 
angle  (J>. 


for  the  density  of  all  points  when  co  >  |x"sin<|)|.  This  density  is 
the  modulation  transfer  function  (MTF)  of  the  ideal  (no  atten¬ 
uation,  perfect  spatial  resolution)  tilted  parallel-beam  measure¬ 
ment.  For  co  <  |x,/sin(|)|  the  MTF  is  zero,  so  those  frequencies 
cannot  be  recovered.  These  frequencies  can  be  represented  as 
missing  conical  regions  whose  volumes  increase  with  tilt  angle 

<t>  (Fig*  4). 

For  co  >  |x"  sin<|>|,  (5)  indicates  that  the  MTF  can  be  inverted 
by  the  scaled  ramp  filter,  |v"cos<j)|.  This  geometrical  derivation 
is  consistent  with  previous  derivations  of  the  backprojection  fil¬ 
ter  for  the  tilted  parallel-beam  geometry  [2]. 


B.  Simulated  Breast  Scan 

In  order  to  compare  OSEM  and  FBP  for  tilted  parallel-beam 
breast  imaging,  a  vertical  axis  of  rotation  (VAOR)  scan  was 
simulated  using  a  section  of  the  digital  anthropomorphic  Pe¬ 
ter  phantom  [11].  This  computer  phantom  consists  of  various 
volumetric  regions  in  the  torso  (Fig.  5).  For  this  study,  only  the 
left  breast  (11  cm  base-to-nipple,  14  cm  base  diameter)  con¬ 
taining  a  uniform  activity  concentration  was  scanned  and  no 
torso  activity  was  present.  The  phantom  was  implemented  on 
a  168x168x168  grid,  voxels  0.35  cm  on  a  side.  The  camera’s 
AOR,  coincided  with  a  line  from  the  base  of  the  breast  through 
the  nipple  (Fig.  6).  The  simulated  VAOR  scan  used  a  22.4  cm  x 
22.4  cm  field-of-view  (FOV)  parallel-hole  gamma  camera,  and 
each  projection  data  set  was  acquired  for  128  angular  views  over 
360°.  Each  projection  was  64x64  pixels  and  pixel  size  was  0.35 
cm  square.  Projection  data  were  simulated  for  0°,  15°,  30°,  and 
45°  tilt  angles.  Expected  values  for  VAOR  projection  data  were 
analytically  computed  based  on  forward  projections  of  activity 
within  the  breast  region  only.  Compton  scattering,  attenuation, 
and  collimator  and  crystal  spatial  resolution  were  not  modeled. 


Area  imaged 


Fig.  5.  Transverse  view  of  digital  anthropomorphic  Peter  phantom.  Squared 
breast  region  indicates  the  subsection  of  phantom  which  was  imaged. 


Camera  head 


Fig.  6.  Illustration  of  simulated  breast  scan.  Breast  represents  subsection  of 
digital  anthropomorphic  Peter  phantom  and  <j>  represents  angular  camera 
tilt. 


An  ensemble  study  was  generated  consisting  of  50  acquisi¬ 
tions  per  tilt  angle  and  based  on  Poisson  variation  about  the 
expected  values.  OSEM  reconstructions  were  performed  with 
the  tilted  geometry  modeled  in  the  photon  detection  probabil¬ 
ities,  and  noise  and  bias  values  were  computed  for  1-50  itera¬ 
tions  (8  subsets).  Reconstructed  OSEM  images  were  32x32x32 
cubic  voxels  each,  0.7  cm  wide.  Compton  scattering,  attenua¬ 
tion,  and  collimator  and  crystal  spatial  resolution  were  not  mod¬ 
eled  within  the  reconstruction.  EBP  reconstruction  was  simi¬ 
larly  performed  with  images  also  32x32x32  cubic  voxels,  0.7 
cm  wide.  EBP  images  were  post-filtered  using  a  3-D  Butter- 
worth  filter.  Butterworth  filter  cut-off  frequencies,  the  frequen¬ 
cies  at  which  the  gain  equals  0.707,  ranged  from  0.1  -  1.0  cy¬ 
cles/pixel  with  frequency  order  equal  to  10. 

Noise  vs.  bias  curves  were  generated  for  both  OSEM  and 
EBP  reconstructions.  Noise  was  computed  as  the  standard  de¬ 
viation  over  the  ensemble  on  a  pixel-by-pixel  basis.  Bias  was 

defined  as  (*,■— Jtf*)2,  where  N  is  the  number  of  voxels, 

x i  is  the  iih  pixel  value  of  the  average  image,  and  is  the  value 

of  the  ith  pixel  in  the  true  phantom  image. 

C.  Defrise  Disk  Phantom  Scan 

Mini-Defrise  phantom  disks  (Data  Spectrum  Corp.,  Hills¬ 
borough,  NC)  were  placed  inside  a  fillable  breast  phantom 
(Model  ECT/FL  BR/A,  Data  Spectrum  Corp.,  Hillsborough, 
NC)  (Fig.  7).  Data  were  acquired  using  one  head  of  a  tiltable 
two  headed  SPECT  system  (Vision  T-22,  Summit  Nuclear, 
Twinsburg,  OH)  and  a  parallel  beam  low  energy  high  resolution 
(LEHR)  collimator  (Hitachi,  Tokyo,  Japan)  (hole  length  =  41.0 


Fig.  7.  Schematic  of  breast  phantom  with  three  mini-Defrise  disks.  Nipple  is 
at  top  and  curvature  at  bottom  corresponds  to  torso  shape. 


mm,  flat  to  flat  hexagonal  hole  size  =  1.8  mm,  septal  thickness 
=  0.18  mm).  The  breast  phantom  contained  three  cold  Defrise- 
phantom  disks  measuring  1.0  cm  in  thickness  and  8.0  cm  in 
diameter  and  with  a  1.0  cm  gap  between  disks.  The  breast  phan¬ 
tom  was  filled  with  water  having  a  22.5  /iCi/ml  concentration  of 
Tc-99m.  The  axis  of  the  breast  phantom  was  then  aligned  par¬ 
allel  to  the  camera’s  AOR  through  the  nipple,  and  the  radius  of 
rotation  (ROR)  was  measured  perpendicular  to  the  center  of  the 
camera  and  to  the  nipple.  Projection  data  were  acquired  for  tilt 
angles  of  0°,  5°,  10°,  and  15°  having  ROR’s  of  12.1  cm,  11.8 
cm,  11.4  cm,  and  10.5  cm  respectively.  The  ROR  decreased 
with  greater  tilt  angle  due  to  greater  access  to  the  breast.  To¬ 
tal  acquisition  time  was  10  seconds  per  angular  step  or  21.33 
minutes  over  360  degrees.  Acquisition  times  were  adjusted  to 
compensate  for  radioactive  decay. 

OSEM  (8  subsets,  3  iterations)  and  FBP  (Butterworth  filtered, 
cut-off  =  0.5  cycles/cm,  order  =10)  reconstructions  were  per¬ 
formed  for  all  tilt  angles,  and  final  images  were  128x128x128 
cubic  voxels,  0.48  cm  on  a  side.  Iteration  and  filter  parameters 
were  determined  from  aforementioned  noise  vs.  bias  study. 

D.  Fillable  Breast  Phantom  with  Lesions 

The  fillable  breast  phantom  (950  ml)  containing  two  spheri¬ 
cal  lesions,  1 .0  cm  diameter  (0.52  ml)  and  0.6  cm  diameter  (0. 1 1 
ml)  (models  ECT/HOL-468/A,  Data  Spectrum  Corp.,  Hillsbor¬ 
ough,  NC),  was  also  imaged  in  this  experiment  (Fig.  8).  Data 
were  acquired  using  a  <  13  cm  x  13  cm  FOV  LumaGEM™ 
gamma  camera  (Gamma  Medica,  Inc.,  Northridge,  CA)  (hole 
length  =  23.6  mm,  flat  to  flat  hexagonal  hole  size  =  1.22  mm, 
septal  thickness  =  0.2  mm)  as  a  component  on  a  dedicated 
Application  Specific  Emission  Tomograph  (ASET)  gantry  [12] 
(Fig.  9).  The  breast  phantom  contained  "mTc-pertechnetate 
(140.6  keV)  activity  and  water  and  lesions :breast  activity  con¬ 
centration  was  5.94:1.  Images  were  acquired  for  128  angular 
views  over  360°,  and  each  projection  was  54x54  pixels,  each 
2.2  mm  square.  Total  acquisition  time  was  20  minutes,  or  9.375 
seconds  per  angular  step.  Projection  data  were  acquired  for  0°, 
10°,  20°,  and  30°  tilt  angles  with  ROR’s  of  7.47  cm,  6.63  cm, 
5.48  cm,  and  5.48  cm  respectively,  measured  from  the  camera 
head  to  the  center-of-rotation  point.  Greater  tilted  angles  allow 
for  the  smaller  ROR’s. 

OSEM  (8  subsets,  3  iterations)  and  FBP  (Butterworth  filtered, 
cut-off  =  0.5  cycles/cm,  order  =  10)  reconstructions  were  per¬ 
formed  for  all  tilt  angles,  and  final  images  were  54x54x54  cubic 
voxels,  0.22  cm  on  a  side.  Iteration  and  filter  parameters  were 
determined  from  aforementioned  noise  vs.  bias  study  to  provide 
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Fig.  8.  Schematic  of  950  ml  breast  phantom  containing  lesions.  Large  lesion 
is  located  at  a  superior  medial  aspect  of  the  breast,  6  cm  from  nipple  and 
3  cm  from  a  normal  to  the  surface  of  the  breast  Smaller  lesion  is  located 
at  an  inferior  medial  aspect  of  the  breast,  6  cm  from  nipple  and  3  cm  from 
a  normal  surface  of  the  breast.  (Top  schematic  shows  relative  height  and 
lesion  displacement;  Bottom  schematic  illustrates  actual  lesion  orientation 
in  the  breast) 


Fig.  9.  Photograph  of  tillable  breast  phantom  suspended  over  the  ASET  system. 


similar  noise  levels  for  both  methods. 

For  both  reconstruction  methods  at  all  degrees  of  tilt,  an  ellip¬ 
tical,  1 1  pixel  region-of-interest  (ROI)  was  taken  for  the  large  le¬ 
sion  and  a  similar  121  pixel  ROI  was  taken  in  the  breast.  Lesion 
contrast  was  calculated  as  ?  and  SNR  as  Pinion- fanast  f 

where  jx  =  mean  and  a  =  standard  deviation  of  selected  ROl’s. 

III.  Results 
A.  Simulated  Breast  Scan 

Figure  10  presents  noise  vs.  bias  plots  for  OSEM  and  FBP 
reconstructions  of  the  simulated  breast  scan.  Figure  1 1  presents 
transaxial,  sagittal,  and  coronal  mid-slices  (with  respect  to  the 
breast  reference  frame)  for  FBP  (no  post-filtering),  and  OSEM 
(3  iterations)  reconstruction  methods  at  0°,  15°,  30°,  and  45°. 
Iteration  cut-offs  were  chosen  to  provide  approximately  equal 
noise  for  both  reconstruction  methods.  FBP  results  indicate  a 
greater  increase  in  backprojection  artifacts  with  greater  angu¬ 
lar  tilt.  OSEM  reconstructions  were  better  able  to  maintain  a 
clear,  defined  border  around  the  breast,  but  both  methods  dis¬ 
played  some  degree  of  axial  stretching  of  the  breast  as  well  as 
inaccurate  activity  distribution  with  greater  tilt.  Stretching  ef¬ 


OSEM 


FBP 


Fig.  10.  Noise  vs.  bias  plots  for  OSEM  (top)  and  FBP(bottom)  reconstructions 
of  simulated  breast  scan  at  0°,  15°,  30°,  and  45°  tilt.  OSEM  data  points 
coirespond  to  iterations  1  through  50  in  increments  of  one  iteration  from 
left  to  right.  FBP  data  points  indicate  Butterworth  filter  cutoffs  (1.0  to 
0.1  cycles/pixel  in  increments  of  0.1  cycles/pixel  from  right  to  left,  where 
furthest  right  point  indicates  no  post-filtering).  OSEM  yielded  the  lower 
bias  for  similar  noise  values. 


fects  were  not  found  transaxially,  however,  due  to  complete  fre¬ 
quency  sampling  of  the  t  =  0  plane  (Fig.  3).  This  is  consistent 
with  tilted  cone-beam  reconstruction  [13],  [14],  where  distor¬ 
tions  were  most  pronounced  along  the  direction  parallel  to  the 
rotation  axis  but  not  along  the  transverse  direction.  Ideally,  the 
breast  would  be  fully  angularly  sampled  to  satisfy  the  Orlov 
condition,  and  hence  reduce  currently  observed  elongation  and 
object  distortion  artifacts  [1],  [4]. 

FBP  noise  vs.  bias  plots  indicate  greater  bias  with  tilt  angle 
compared  to  OSEM  reconstructions  (Fig.  10).  A  50%  increase 
in  bias  was  noted  between  FBP  0°  and  15°  recontructions,  while 
only  a  15%  increase  in  bias  was  noted  between  OSEM  0°  and 
15°  reconstructions,  compared  at  similar  noise  levels.  Butter- 
worth  filtering  tended  to  offer  the  best  compromise  between 
noise  and  bias  for  FBP  reconstructions,  specifically  at  a  0.3  cy¬ 
cles/pixel  cutoff.  Three  iterations  displayed  the  best  noise/bias 
trade-off  for  OSEM  0°,  15°,  and  30°  data  and  four  iterations 
were  found  optimal  for  45°  data.  OSEM  reconstructions  dis¬ 
played  the  lower  overall  bias. 

B.  Defrise  Disk  Phantom  Scan 

Figure  12  displays  reconstructions  of  mini-Defrise  phantom 
disks  placed  inside  a  tillable  breast  phantom.  The  resulting  pro¬ 
files  indicate  a  bowing  effect  slightly  visible  at  10°  tilt,  and 
more  clearly  visible  at  15°  tilt.  OSEM  Defrise  disk  contrast 
was  significantly  better  at  greater  tilt  angles  compared  to  FBP 
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Fig.  11.  Transaxial,  coronal,  and  sagittal  (with  respect  to  the  breast  reference  frame)  OSEM  and  FBP  reconstruction  slices  for  simulated  breast  scan.  Slices  were 
taken  for  0°,  15°,  30°,  and  45°  tilt  and  corresponding  truth  phantom  slices  located  at  bottom. 


TABLE  I 

Contrast  and  SNR  values  of  large  lesion  for  OSEM  and  FBP 

RECONSTRUCTIONS. 

0°  Tilt 

SNR  Contrast  ROR  (cm) 


OSEM 


0° 

10° 

20° 

30° 

9.98 

13.43 

12.23 

12.85 

2.65 

3.84 

2.82 

3.57 

7.47 

6.43 

5.48 

5.48 

5°  Tilt 

FBP 

10°  Tilt 

0° 

5.57 

2.26 

7.47 

10° 

6.39 

2.13 

6.43 

20° 

6.59 

2.75 

5.48 

w 

e 

0 

9.10 

3.25 

5.48 

15°  Tilt 

OSEM 


FBP 


results.  OSEM  reconstruction  sets  displayed  a  one  pixel  (0.48 
cm,  3.8%)  increase  in  breast  length  at  15°  tilt,  measured  from 
the  nipple  to  the  anterior  wall.  FBP  profiles  indicated  a  two 
pixel  (0.98  cm,  7.6%)  increase  in  breast  length.  All  profiles 
drawn  axially  across  the  right-center  of  the  breast  indicate  an 
increase  in  counts  toward  the  nipple  due  to  attenuating  effects 
of  the  Defrise  disk  phantom  within  the  breast. 


C.  Fillable  Breast  Phantom  with  Lesions 

Table  I  presents  contrast  and  SNR  values  for  reconstruc¬ 
tions  of  the  fillable  breast  phantom  containing  lesions  (Fig.  13). 
OSEM  reconstructions  showed  greater  SNR  and  contrast  values 
at  all  tilt  angles.  Because  camera  tilt  allows  greater  access  to  the 
breast  and  in  turn  decreases  the  effective  ROR,  overall  contrast 
and  SNR  values  generally  increased  with  tilt  angle,  confirming 
the  effectiveness  which  can  be  achieved  with  TH-SPECT  and 
dedicated  compact  gamma  cameras. 


Fig.  12.  Sagittal  view  OSEM  and  FBP  reconstructions  of  mini-Defrise  phantom 
combined  with  tillable  breast  phantom  and  corresponding  profiles.  Profile 
markers  are  drawn  axially  across  left-center  of  breast  phantom. 
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Fig.  13.  Sagittal  view  OSEM  and  FBP  reconstructions  of  anthropomorphic 
torso  phantom  combined  with  fillable  breast  phantom  containing  lesions. 
Profile  markers  are  drawn  axially  across  left-center  of  breast  phantom. 


IV.  Discussion  and  Conclusion 


riety  of  different  3-D  imaging  systems  with  limited  view  angles 
[2],  [15],  [16],  [17].  If  frequency  components  within  the  miss¬ 
ing  cone  and  especially  along  the  x  axis  are  assumed  to  be  zero, 
the  reconstructed  mean  densities  of  all  planes  perpendicular  to  z 
will  be  equal.  This  has  the  effect  of  producing  unwanted  activity 
in  areas  where  activity  should  be  zero  or  near-zero.  These  activ¬ 
ity  distortions  are  observed  in  FBP  techniques  as  streak  artifacts 
(Fig.  11)  along  the  axis  of  rotation.  Due  to  complete  sampling 
of  the  x  =  0  plane,  transverse  slices  do  not  exhibit  streak  arti¬ 
facts,  although  in-plane  distortions  are  present. 

The  quality  of  TH-SPECT  reconstructions  clearly  degrades 
with  tilt  angle,  specifically  due  to  the  volumetric  increase  in  null 
space.  The  fraction  of  null  space  to  total  frequency  volume  is 
given  as 


Null  space 
Total  volume 


=  sin(J). 


(6) 


With  30°  tilt,  only  13.4%  of  frequency  space  accounts  for 
missing  data.  With  60°  tilt,  however,  half  of  the  frequency  data 
is  missing.  This  non-linear  increase  in  unsampled  volume  indi¬ 
cates  tilt  angles  over  30°  may  not  be  useful  for  breast  SPECT. 

The  conjugate  views  for  non-tilted  projections  exist,  and 
therefore  0  need  only  be  sampled  from  0  to  n.  For  tilted  projec¬ 
tions,  however,  the  conjugate  view  does  not  exist  and  0  should 
be  sampled  from  0  to  271.  This  is  illustrated  in  Figure  3  (bottom) 
as  rotation  from  0°  to  180°  degrees  does  not  even  fully  sample 
the  frequency  space  region  outside  the  circle  of  radius  x"  sin<|>. 

Herein  we  have  investigated  tilted  orbits  in  which  the  camera 
center  moves  along  a  circle,  elucidating  the  nature  of  the  incom¬ 
plete  sampling  provided  by  these  orbits  and  the  effectiveness  of 
FBP  and  OSEM  in  utilizing  this  incomplete  data.  In  other  work, 
we  proposed  circle-plus-arc  [1],  [4]  and  more  elaborate  orbits 
that  do  provide  complete  sampling  of  an  isolated  breast,  and  we 
have  developed  a  SPECT  system  capable  of  implementing  such 
orbits  [12]. 
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Abstract 

As  a  test  bed  for  dedicated  breast  SPECT  cameras 
that  are  under  development,  a  general  purpose  SPECT 
system  with  two  tiltable  heads  was  used  to  image  Tillable 
breast  and  torso  phantoms  containing  multiple  lesions. 
Breast,  liver,  and  myocardial  activity  were  included  in  order 
to  simulate  direct  contamination  and  Compton  scattering 
expected  in  clinical  scans.  The  TH-SPECT  data  were 
reconstructed  using  an  OS-EM  algorithm  which  accounted 
for  the  tilted  geometry.  High  count  planar  images  were 
acquired  for  comparison  with  tilted-head  SPECT  (TH- 
SPECT).  In  order  to  characterize  axial  blurring  effects 
inherent  with  TH-SPECT  reconstructions,  two  cylindrical 
disk  Defrise  phantoms,  one  large  Defrise  phantom  and  one 
mini-Defrise  phantom  placed  inside  the  Tillable  breast 
phantom,  were  imaged  at  various  tilt  angles.  Results 
indicate  an  increase  in  axial  blurring  with  greater  tilt  angle. 
Reconstructions  of  the  combined  Tillable  breast  and  torso 
phantoms  containing  two  1.15  ml  lesions,  one  centered 
axially  and  one  proximal  to  the  anterior  chest  wall  within 
the  breast,  were  most  clearly  visible  in  the  30°  reconstructed 
TH-SPECT  images,  providing  lesion  SNR  and  contrast 
improvements  of  nearly  three  times  compared  to  the  high- 
count  planar  images 

I.  INTRODUCTION 

For  breast  cancer  imaging  using  conventional  Single 
Photon  Emission  Computed  Tomography  (SPECT)  systems, 
the  gamma  camera’s  axis  of  rotation  is  parallel  to  the  patient 
table  -  a  horizontal  axis  of  rotation  (HAOR).  This  geometry 
decreases  breast  signals  due  to  torso  attenuation  and,  more 
importantly,  requires  a  large  radius  of  rotation  which  yields 
poor  results  for  SPECT  breast  imaging  [1]. 

Previous  work  indicates  breast  cancer  imaging  can 
be  improved  by  using  a  vertical  axis  of  rotation  (VAOR)  [1], 
[2].  In  this  work,  we  expand  upon  previous  VAOR  breast 
imaging  by  using  tiltable  gamma  cameras.  By  tilting  the 
camera  off  axis  with  respect  to  the  nipple-to-anterior  chest 
axis  through  the  breast,  improved  imaging  of  lesions  near 
the  chest  wall  and  a  reduced  radius  of  rotation  (ROR)  can  be 
achieved.  This  paper  investigates  tilted  VAOR  breast 
imaging  (Fig.  1)  by  simulating  a  dedicated  VAOR  system 
using  a  HAOR  SPECT  system  with  tiltable  gamma  cameras. 

One  major  concern  with  Tiltable-Head  SPECT 
(TH-SPECT)  imaging  is  insufficient  angular  sampling.  The 
completeness  of  reconstructions  using  TH-SPECT  is  bound 
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by  Orlov  sphere  conditions  [3].  In  order  to  satisfy  these 
conditions,  the  image  space  must  be  sampled  over  a  great 
circle  or  equatorial  orbit.  By  tilting  the  camera  heads  off 
axis,  this  primary  requirement  is  not  met  (Fig.  2),  and  thus 
there  is  the  potential  for  artifacts  in  the  reconstructed  images. 
This  is  also  an  issue,  in  general,  with  slant-hole  collimators 
and  rotating  slant-hole  collimators  [4]. 

Another  concern  with  TH-SPECT  imaging  is 
background  interference  from  myocardium,  liver,  and  body 
tissue  that  is  expected  clinically  with  the  use  of  Tc-99m 
labeled  tracers.  By  tilting  the  camera  heads  off  axis,  heart 


Figure  1 :  Schematic  of  dedicated  breast 
VAOR  system. 


Figure  2:  Sphere  showing  an  equatorial  camera 
orbit  and  an  orbit  where  the  camera  head  is  tilted. 
The  camera  is  normal  and  tangentral  to  the  sphere 
at  all  projection  views. 


and  liver  activity  are  introduced  into  the  camera’s  field  of 
view  (FOV)  at  some  projections,  potentially  contaminating 
TH-SPECT  reconstructions. 

The  following  measured  results  indicate,  however, 
that  myocardium,  liver,  and  tissue  background  interference 
do  not  markedly  affect  breast  lesion  contrast  or  signal-to- 
noise  ratio  (SNR)  values.  Furthermore,  tilted  head 
acquisitions  can  be  successfully  reconstructed,  and  although 
axial  distortions  are  present  and  significantly  increase  at 
greater  tilt  angles,  excellent  SNR  and  contrast  values  are 
obtained  for  lesions,  including  those  near  the  chest  wall. 

II.  MATERIALS  AND  METHODS 

All  data  were  acquired  using  a  tiltable  two-headed 
SPECT  system  (Vision  T-22,  Summit  Nuclear,  Twinsburg, 
OH)  (Fig.  5).  Only  one  head  was  used  to  obtain  data  during 
this  study.  A  parallel  beam  low-energy  high-resolution 
(LEHR)  collimator  (Hitachi)  was  used  in  all  experiments. 
Collimator  hole  length  was  41.0  mm,  hexagonal  hole  size 
was  1.8  mm,  flat-to-flat  and  septal  thickness  was  0.18  mm. 
All  images  were  reconstructed  by  OS-EM  (5  iterations,  8 
subsets),  with  the  tilted  geometry  modeled  in  the  photon 
detection  probabilities.  Compton  scattering,  attenuation,  and 
collimator  and  crystal  spatial  resolution  were  not  modeled. 

A.  Defrise  Disk  Phantom 

A  Defrise  disk  phantom  (Fig.  3)  (Data  Spectrum 
Corp.,  Hillsborough,  NC)  was  utilized  in  the  first  experiment. 
This  phantom  consisted  of  a  cylinder  containing  six  acrylic 
disks  placed  1.5  cm  apart.  The  acrylic  disks  had  a  thickness 
of  1.2  cm  and  a  diameter  of  18.9  cm.  Five  emission  disks 
were  then  formed  by  adding  16.6  mCi  of  Tc-99m  and  water 
to  the  cylinder,  each  hot  disk  having  a  uniform  concentration 
of  4.1  pCi/ml.  The  phantom  was  then  placed  parallel  to  the 
camera’s  axis  of  rotation  (AOR)  with  a  radius  of  rotation 
(ROR)  of  11.5  cm,  measured  from  the  center  of  the  camera 
to  the  Defirse  disk  phantom’s  long  axis,  perpendicular  to  the 
camera’s  AOR  (Fig.  3).  Projection  data  were  taken  for  0°, 
15°,  30°,  45°,  and  60°  tilt  angles.  Each  projection  data  set 
was  acquired  for  128  angular  views  over  360°,  and  each 
projection  was  128x128  pixels  and  pixel  size  was  0.48  cm 
square.  Total  acquisition  time  was  7  seconds  per  angular 
step  for  a  total  time  of  15  minutes.  Subsequent  acquisition 
times  were  increased  to  compensate  for  radioactive  decay. 

B.  Breast  Phantom  with  Mini-Disk  Phantom 

In  order  to  characterize  axial  blurring  within  the 
breast  geometry,  custom  mini-Defrise  phantom  disks  (Data 
Spectrum  Corp.,  Hillsborough,  NC)  were  placed  inside  a 
custom  Tillable  breast  phantom  (Data  Spectrum  Corp., 
Hillsborough,  NC)  (Fig.  4).  By  filling  the  breast  phantom 
with  water  having  a  22.5  pCi/mL  concentration  of  Tc-99m, 
three  cold  disks  measuring  1.0  cm  in  thickness  and  8.0  cm  in 
diameter  and  two  hot  disks  measuring  1.0  cm  and  8.0  cm  in 
diameter  were  formed.  The  axis  of  the  breast  phantom  was 
then  placed  parallel  to  the  camera’s  AOR  and  the  radius  of 
rotation  was  measured  from  the  nipple  to  the  center  of  the 
camera,  perpendicular  to  the  camera’s  AOR.  Because  axial 
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Figure  3:  Schematic  of  the  Defrise  phantom 
with  camera  head  at  0°,  15°,  30°,  45°,  and  60° 
tilt  angles.  Curved  arrow  illustrates  camera 
orbit  out  of  page  about  the  AOR. 


Figure  4:  Photograph  of  the  custom  breast  phantom  with 
three  mini-Defrise  disks  8.0  cm  in  diameter  and  1.0  cm 
thick,  separated  by  a  1.0  cm  gap  and  aligned  along  the 
breast  axis.  Nipple  is  at  right,  and  curvature  at  left 
corresponds  with  torso  shape. 


blurring  becomes  extremely  noticeable  at  larger  tilt  angles 
(>  15°),  projection  data  were  taken  for  angles  of  0°,  5°,  10°, 
and  15°  having  ROR’s  of  12.1  cm,  11.8  cm,  11.6  cm,  and 
10.9  cm  respectively.  The  ROR  was  decreased  with  greater 
tilt  angle  due  to  greater  access  to  the  breast.  Potentially 
smaller  ROR’s  would  be  possible  with  dedicated  small  area 
breast  imaging  systems.  Total  acquisition  time  was  10 
seconds  per  angular  step  or  21.33  minutes.  Subsequent 
acquisition  times  were  increased  to  compensate  for 
radioactive  decay. 

C.  Anthropomorphic  Torso  Phantom  and  Breast 
Phantom  with  Lesions 

An  anthropomorphic  torso  phantom  (Model 
ECT/TOR/P,  Data  Spectrum  Corp.,  Hillsborough,  NC) 


Figure  5:  Photograph  of  anthropomorphic  torso  and 
breast  phantoms  with  dual  tilted  heads  on  the  SPECT 
camera. 

combined  with  the  tillable  breast  phantom  was  imaged  in  this 
experiment  (Fig.  5)  to  simulate  the  expected  emission 
contamination  from  the  myocardium  and  liver  to  the  specific 
breast  imaging  task.  The  craniocaudal  axis  of  the  torso 
phantom  was  perpendicular  to  the  camera’s  axis  of  rotation, 
thus  necessitating  the  tilted  camera  heads  and  allowing  for 
close  proximity  imaging  of  the  breast  phantom.  Two  1.3  cm 
diameter  lesions  (1.15  ml)  with  Tc-99m  concentration  of 
11.1  p.Ci/ml  were  suspended  inside  the  left  breast.  One 
lesion  was  placed  axially  in  the  center  of  the  breast,  5.1  cm 
from  the  nipple.  The  second  lesion  was  located  2.2  cm 
anterior  to  the  chest  wail  and  positioned  3.0  cm  towards  the 
medial-superior  side  of  the  breast.  A  third  0.9  cm  diameter 
lesion  (0.38  ml)  with  a  concentration  of  11.1  pCi/ml  was 
placed  outside  the  breast  on  the  chest  wall  to  represent 
axillary  tumor  growth.  The  breast  phantom  had  a  uniform 
concentration  of  1.5  pCi/ml.  In  order  simulate  contamination 
effects  due  to  direct  superposition  of  myocardial  activity  on 
the  breast  at  some  projection  views  and  also  Compton 
scattering  expected  in  clinical  scans,  the  tumor  and  breast 
were  imaged  with  additional  liver  and  myocardium  activity. 
The  lesion:liver:myocardium:breast-background  ratio  was 
7.4:12.5:12.5:1,  giving  the  liver  and  heart  a  Tc-99m 
concentration  of  18.7  pCi/ml. 

Projection  data  were  acquired  for  0°,  15°,  30°,  and 
45°  tilt  angles  with  ROR’s  of  18.1  cm,  10.5  cm,  9.4  cm,  and 
5.9  cm  respectively,  measured  as  previously  described.  The 
images  were  acquired  for  128  angular  views  over  360°,  and 
each  projection  was  128x128  pixels  and  pixel  size  was  0.48 
cm  square.  Total  acquisition  time  was  32  minutes,  or  15 
seconds  per  angular  step  and  subsequent  acquisition  times 
were  increased  to  compensate  for  radioactive  decay. 

Because  planar  scintimammography  is  gaining 
acceptance  as  an  effective  way  to  identify  breast  lesions 
greater  than  1  cm  in  diameter  [5],  a  planar 
scintimammography  image  (craniocaudal  view)  was  made 
for  comparison.  Breast,  liver,  and  myocardium  background 
radiation  were  present,  with  the  same  activity  concentration 
as  for  TH-SPECT,  and  the  total  acquisition  time  was  5 
minutes. 


30  °  45  ° 
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Figure  6:  Single  slice  sagittal  views  of  OS-EM 

reconstructed  Defrise  phantom  images  at  indicated  tilt 
angles. 

III.  RESULTS 

A .  Degradative  Angular  Sampling  Effects 

Single  Slice  sagittal  views  of  the  reconstructed 
Defrise  phantom  images  at  0°,  15°,  30°,  45°,  and  60°  are 
presented  in  Figure  6.  For  the  reconstructed  sagittal  slices, 
each  disk  should  be  regular  "bar"  shape  (Fig.  6.,  left,  top), 
but  as  the  tilt  angle  is  increased  there  is  a  clear  bowing  effect 
in  the  axial  direction  which  corroborates  the  results  in  [6]. 
The  disks  are  almost  unrecognizable  in  the  30°-60°  tilts. 

Figure  7  compares  the  results  of  the  mini-Defrise 
phantom  inserted  into  the  tillable  breast  phantom.  Results 
from  the  larger  Defrise  phantom  indicated  extreme  axial 
blurring  for  tilt  angles  greater  than  15°.  This  experiment 
investigated  smaller  tilt  angles  to  evaluate  the  effects  of 
incomplete  sampling  with  small  tilt  angles,  and  also  to 
attempt  to  minimize  the  blurring  effect  but  still  maximize 
chest  wall  imaging  without  additional  orbits  to  complete  the 
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Figure  7:  Single  slice  sagittal  view  OS-EM 
reconstructions  of  mini-Defrise  phantom  combined  with 
fillable  breast  phantom  at  indicated  angles.  Profile 
markers  are  drawn  axially  across  left-center  of  breast 
phantom. 


Figure  8:  Profiles  drawn  across  single  slice  sagittal  views  of 
mini-Defrise  phantom  (Fig.  7).  Left  corresponds  to  bottom 
portion  of  the  images  shown  in  Figure  7.  Camera  head  was  tilted 
away  from  chest  wall.  Note  increase  in  breast’s  axial  length 
with  greater  tilt  angle. 
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Figure  9:  Sagittal  view  reconstructions  of  fillable  breast  phantom 
attached  to  anthropomorphic  torso  phantom  containing  two  hot 
lesions  (one  lesion  axially  centered,  one  lesion  anterior  to  the 
chest  wall  and  positioned  towards  the  medial-superior  side  of  the 
breast)  at  0°,  15°,  30°,  and  45°  tilt  angles  as  well  as  conventional 
high-count  planar  images.  Profiles  for  each  lesion  are  also 
shown.  The  cold  spot  directly  behind  the  central  lesion  is  the 
acrylic  support  used  to  hold  lesion  in  place  and  the  cold  gap 
between  the  breast  and  torso  represents  phantom  chest  wall. 


Table  1 


SNR  and  contrast  values  for  center  lesion  and  offset  lesion  at 
0°,  15°,  30°,  and  45°  tilt  angles  as  well  as  a  conventional  high- 
count  planar  image  (craniocaudal  view). 


SNR 

[(^1  lesion  —  JJ-bieest  )/  CTbrcast] 

Contrast 

[(glesion  ”  |^bre»s»)/  ^breast] 

ROR 

(cm) 

Center  lesion 

Planar 

2.10 

0.20 

18.1 

0° 

2.09 

0.25 

18.1 

15° 

3.56 

0.49 

10.5 

30° 

4.89 

0.59 

9.4 

45° 

5.15 

0.59 

5.9 

Offset  lesion 

Planar 

2.04 

0.20 

18.1 

0° 

2.13 

0.24 

18.1 

15° 

3.74 

0.51 

10.5 

30° 

6.32 

0.63 

9.4 

45° 

5.16 

0.49 

5.9 

data.  A  bowing  effect  similar  to  that  seen  in  Figure  6  is 
clear  in  these  reconstructions,  although  not  as  pronounced. 

Profiles  drawn  axially  across  the  left-center  of  the 
breast  (Fig.  8)  indicate  an  increase  in  counts  towards  the 
nipple  due  to  the  attenuating  effects  of  the  defrise  phantom 
within  the  breast.  The  profiles  also  show  a  definite  axial 
"stretching"  effect,  as  the  15°  image  is  roughly  0.5  cm  longer 
than  the  0°  image.  This  result  indicates  that  visualization  of 
axially  juxtaposed  tumors  or  shell-like  tumors  may  be 
difficult  to  discern. 

B.  Lesion  Imaging  with  Various  Tilt  Angles 

Reconstructed  sagittal  views  of  the  Tillable  breast 
phantom  containing  two  hot  lesions  combined  with  the 
anthropomorphic  torso  phantom  are  presented  in  Figure  9 
with  corresponding  profiles  at  the  indicated  locations. 
Background  heart,  liver  and  torso  contamination  were 
present,  and  all  profiles  were  drawn  across  the  center  of  the 
lesions.  Table  1  compares  the  SNR  and  contrast  values  for 
each  lesion  and  tilt  angle.  Figure  9  clearly  shows  an  increase 
in  both  contrast  and  SNR  as  tilt  angle  is  increased  from  0°  to 
30°.  Tilt  angles  beyond  30°  do  not  significantly  improve 
either  contrast  or  SNR  values. 

The  third  lesion,  attached  to  the  anthropomorphic 
torso  phantom  anterior  to  the  chest  wall  and  superior  to  the 
left  breast,  was  also  reconstructed.  Although  not  shown,  the 
lesion  is  clearly  visible  in  all  reconstructions. 


IV.  DISCUSSION  AND  CONCLUSIONS 

Because  of  insufficient  sampling,  axial  blurring  is 
evident  in  all  TH-SPECT  reconstructions.  As  shown  with 
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Figure  10:  Transverse  view  of  phantom  illustrating 
TH-SPECT  acquisition  as  well  as  an  additional  HAOR 
orbit  to  gather  insufficiently  sampled  data. 

the  larger  Defrise  phantom,  large  tilt  angles  (45°  and  60°) 
yield  extremely  distorted  images.  Because  the  ROR  can  be 
decreased  as  tilt  angle  is  increased,  it  is  important  to  find  an 
acceptable  combination  between  the  two,  specifically  one 
which  minimizes  axial  blurring  effects  while  still 
maximizing  SNR  and  contrast  values.  Trade-offs  between 
higher  resolution  parallel  hole  collimators  versus  higher 
sensitivity  collimators  should  be  considered.  The  use  of 
non-parallel  hole  collimators  should  also  be  considered  in 
this  application  since  these  geometries  could  potentially  help 
complete  undersampled  regions  of  the  Orlov  sphere  [7]. 

As  Figure  9  shows,  larger  tilt  angles  can  "stretch" 
the  breast  axially.  This  stretching  phenomena  also  affects  the 
lesion,  and  can  potentially  blur  together  axially  offset  lesions 
or  distort  shell-like  lesions.  This  phenomena  has  very  little 
effect,  however,  on  the  15°  degree  reconstruction. 
Compared  to  the  0°  reconstruction  which  has  no  distortions, 
SNR  and  contrast  values  between  the  two  are  significantly 
different,  providing  lesion  SNR  and  contrast  improvements 
of  up  to  two  times  for  the  center  lesion  (Table  1).  Although 
30°  and  45°  tilt  reconstructions  have  even  better  SNR  and 
contrast  values,  distortions  of  the  breast  are  much  more 
evident,  compared  to  the  smaller  tilt  angles. 

As  Table  1  indicates,  effective  radius  of  rotation 
significantly  decreases  as  tilt  angle  is  increased.  Although 
the  45°  radius  of  rotation  was  almost  half  that  of  the  30°, 
extreme  insufficient  sampling  as  well  as  attenuation 
prevented  further  contrast  and  SNR  improvements.  Dramatic 
decreases  in  ROR  due  to  increasing  tilt  angle  could  not, 
however,  be  achieved  using  a  small-area  dedicated  breast 
SPECT  system,  although  results  indicate  SNR  and  contrast 
values  should  improve  with  any  reduction  of  ROR.  In 
addition  to  viewing  lesions  near  the  torso-breast  interface, 
TH-SPECT  reconstructions  may  considerably  improve  SNR 
and  contrast  values  even  for  lesions  within  the  breast, 
presumably  by  allowing  for  a  reduced  ROR.  Localization  of 
the  lesion,  however,  may  be  difficult  with  the  present 
artifacts. 

The  blurring  and  distortion  artifacts  become  more 
apparent  with  greater  axial  non-uniformity  in  the  breast,  as 


shown  in  Figure  7  with  the  mini-Defrise  phantom.  The  15° 
tilt  reconstruction  of  the  mini-Defrise  phantom  inside  the 
Tillable  breast  phantom  highlights  a  definite  bowing  effect. 
These  results  indicate  that  artifacts  will  be  more  present 
when  objects  are  axially  aligned  and  imaged  with  a  TH- 
SPECT  technique.  This  bowing  effect  is  not  visualized  in 
Figure  9  reconstructions  simply  due  to  the  uniformity  of 
breast  tissue  and  spherical  lesions.  Non-uniform  breast 
lesions,  such  as  distorted  or  shell-like  tumors,  may  result  in 
further  image  degradation  even  at  small  angular  tilts. 
Although  dedicated  breast  SPECT  systems  under 
development  may  also  suffer  from  the  artifacts  illustrated  in 
this  study,  those  systems  may  further  provide  SNR  and 
contrast  improvements  due  to  reduced  ROR’s. 

When  the  reconstruction  was  limited  to  the  breast 
only,  essentially  eliminating  background  interference  from 
the  FOV,  a  cupping  artifact  [8]  was  present.  This  cupping 
artifact  is  an  artificial  increase  in  activity  near  the  edge  of  the 
FOV.  By  reconstructing  the  entire  image  space  (heart  and 
liver  interference  included),  this  artifact  was  eliminated  [8]. 
It  remains  to  be  seen  if  small-area  dedicated  breast  SPECT 
systems  will  suffer  from  these  limitations  and  if 
reconstructed  images  from  a  small-area  camera  contain  the 
cupping  artifact. 

Possible  modifications  can  be  made  to  the  TH- 
SPECT  geometry  to  account  for  insufficient  sampling  and  to 
reduce,  if  not  eliminate,  axial  blurring  artifacts.  Because 
Orlov  Sphere  requirements  are  not  met  with  any  degree  of 
camera  tilt  [3],  an  additional  orbit  is  required  to  satisfy  these 
conditions  for  fully  quantitative  SPECT  imaging  [9].  Figure 
10  introduces  one  such  additional  orbit.  If  TH-SPECT 
VAOR  data  is  collected  and  combined  with  an  incomplete 
HAOR  orbit  [9],  Orlov  Sphere  requirements  should  be  met 
and  blurring  due  to  insufficient  axial  sampling  should  be 
eliminated.  This  short  arc-orbit  (ranging  from  10°  to  90° 
depending  upon  tilt  angle)  needs  only  to  gather  enough  data 
to  fully  sample  the  breast,  as  accurate  heart  and  liver 
reconstruction  may  not  be  necessary. 
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